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ABSTRACT 

Lindblad resonances have been suggested as an important mechanism for angular mo- 
mentum transport and heating in discs in binary black hole systems. We present the 
basic equations for the torque and heating rate for relativistic thin discs subjected 
to a perturbation. The Lindblad resonance torque is written explicitly in terms of 
metric perturbations for an equatorial disc in a general axisymmetric, time-stationary 
spacetime with a plane of symmetry. We show that the resulting torque formula is 
gauge-invariant. Computations for the Schwarzschild and Kerr spacetimes are pre- 
sented in the companion paper. 

Key words: accretion, accretion discs - relativistic processes - black hole physics. 



1 INTRODUCTION 

The past several years have seen a surge in interest re- 
lated to the electromagnetic signatures of merging black 
holes. Such a signature would have to come not from 
the black holes themselves, but from the gas that sur- 
rounds them. Heating of this gas and consequent emis- 
sion of electromagnetic radiation h as been discussed both 
in the conte xt of the inspiral phas e (|Chang et al.|[2uTol ). the 
coalescence (|Kocsis fc Lo eb 2008), and in the post-merger 
phase as the mass loss and kick of the final black hole 
modi f y the orbits of the gas particles |B odc & P hinneyl 
20071; IShields fc Bonninel 120081; ISchnittman fc Krolikl |200S| ; 



Anderson et al.ll2010l ; iRossi et alj|2010h . 

It is often suggested that torques arising from 
Lindblad resonance^J play a key role in redistributing 
gas in the inspiral phase |Armitage fc Nataraian 20021; 



Milosavlievic fc Phinnevl [2005T MacFadven fc Milosavlievid 
20081 : IChang et al.ll2010l ) and controlling the surface density 
profile and heating rate of the gas disc. These torques act by 
exciting density perturbations at the location of either inner 
or outer Lindblad resonances (ILRs or OLRs), at which the 
synodic period (i.e. the time between successive passages of 
the secondary black hole and a disc particle) is an integer 
multiple of the period of radial epicyclic oscillations in the 
disc. In some scenarios, the resonant torques operate in the 
nonrelativistic Newtonian regime, which has a long history 



1 Other resonances may also be relevant, e.g. it has been sug- 
gested that there could be matter at the L4 and L5 Lagrange 
points of binary black holes llSchnittmaitolOh , but they require 
a fundamentally different treatment and will not be investigated 
here. 



of study in the context of g alactic discs, planetary rings 



and circumstellar discs (e.g. Lvndcn-Bcll & Kalnais! 19721: 



iGoldreich fc Tremaind [l978l. Il979. 1980l; iLin fc Papaloizoul 



1 19791) . However, in other s - particularly the cases of inner 
discs l|Chang et al.|[2uTol ) - Lindblad resonant torques are 
used all the way in to radii of a fewxlOM. In these cases, 
it is desirable to revisit the Lindblad resonances in a fully 
relativistic context. This is especially true since pericentre 
precession introduces an additional ILR (the m = 1 or 0:1 
ILR) that has no analogue in the Newtonian-Keplerian prob- 
lem. The principal purpose of this paper and its companion 
is to provide a relativistic treatment of the Lindblad torques, 
including computation of the torque formula in black hole 
spacetimes (Schwarzschild or Kerr), in the extreme mass ra- 
tio limit. 

This paper and its companion are not concerned with a 
full analysis of any one scenario for the generation of an elec- 
tromagnetic counterpart to a black hole mer ger, although 
they are most relevant to the proposal of IChang et al.l 
(2010). Rather, our motivation is to establish the relativistic 
Lindblad torque formula so that it can be used to establish 
the role (or lack thereof) of Lindblad torques in future work. 
In this paper ("Paper I"), we develop the general formalism 
for Lindblad torques in thin discs orbiting in the equatorial 
planes of axisymmetric, time-independent spacetimes with 
a plane of symmetry, and with weak perturbations of gen- 
eral form respecting the equatorial reflection symmetry. This 
covers the case of a binary Schwarzschild black hole with an 
extreme mass ratio (q = M2/M1 <C 1) and a gas disc or- 
biting in the same plane. It also covers the Kerr case if the 
primary hole's spin is aligned with the orbital angular mo- 
mentum of the binary and disc (which may or may not be 
the physical case; here it is a simplifying assumption that we 
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may wish to remove in future work) . We work out the torque 
formula in terms of the background metric and its perturba- 
tion h a p and establish generic features such as gauge invari- 
ance. The companion ("Paper II") focuses on the specific 
cases of interest - the Schwarzschild and Kerr metrics with 
a small perturber - and describes the numerical evaluation 
of the resonant torque. 

Our analysis considers the case of geometrically thin 
discs. The alternative - a geometrically thick disc, such as 
that in an advection-dominated accretion flow (ADAF) - 
cannot be treated by the methods described here. A more 
appropriate model for an extreme mass ratio binary where 
the secondary orbits within a thick disc was considered by 
iNaravanl (|2000h . We note however that one conceivable way 
to produce such a thi ck disc, even at initi ally high accretion 
rates as considered bv lChang et al.l (|2010T ). would be for res- 
onant heating to destroy the thin disc solution and result in 
a radiatively inefficient inner disc. Assessment of this possi- 
bility requires us to be able to quantitatively compute the 
resonant torques. 

We evaluate the torque here by assuming a particle disc, 
since previous works on Lindblad resonances have found that 
the specific dissipation mechanism (e.g. viscosity or propa- 
gation of spiral density waves as occurs in a hydrodynamic 
disc) does not affect the total torque at a resonance so long as 
the excitation of disc modes is localized near the resonance 
and in the linear regime fe.g.lMever-Vernet fc Sicardvll 19871 ; 
iLubow fc Qgilvielll998l ; IOgilviell2007l l. The underlying rea- 
son for this - namely, that the vector eccentricitjQ integrated 
over the resonance in each sector of the disc, which is both 
excited by the external perturbation and acted upon by the 
perturbation to yield the overall torque, is not changed but is 
simply redistributed by short-range interactions among disc 
particles - is generic and we expect it to also hold in the 
relativistic case. We also note that the modes of oscillation 
of relativistic discs have been investigated (e.g. |Perez et alj 
1 19971 ; ISilbergleit et all 1200 ll ; lOrtega-Rodriguez et al.ll2002l ). 
however their excitation by perturbations to the spacetime 
have not yet been treated. 

This paper is organized as follows. Section [2] lays out 
the assumed background spacetime and the motion of test 
particles in it. Section [3] describes the behaviour of particles 
under a general perturbation to their Hamiltonian (gravi- 
tational or otherwise) and the resulting torque on an ini- 
tially axisymmetric disc. Section U re-expresses this torque 
in terms of the metric perturbation and demonstrates gauge 
invariance; it also gives a useful alternative expression for 
the torque in terms of the power delivered to a test particle 
on a slightly eccentric orbit. Section [5] shows that our expres- 
sion reduces to the familiar expression for Lindblad torques 
in the familiar Newtonian-Keplerian case, i.e. in the space- 
time of a point mass at radii r/M ^> 1. Section [6] describes 
the disc heating at the resonance, and Section [7] concludes. 

We use relativistic units where G = c = 1. 



2 The vector eccentricity is the eccentricity weighted by the di- 
rection of the pericentre, or ee ln7 in Keplerian elements. In a non- 
Keplerian potential the longitude of pericentre precesses, but the 
vector eccentricities of particles at the same epoch may still be 
summed. 



2 BACKGROUND SPACETIME AND 
PARTICLE TRAJECTORIES 

2.1 The spacetime 

We consider the unperturbed problem of a disc orbiting in 
the equatorial plane of a black hol e. In the equatorial pla ne, 
the metric may be written as (e.g. iPage fc Thornelll974l ) 



-e at + e 



2i/' 



(d<j> - udt) 2 + e 2,i dr 2 + dz 2 



(1) 



where v, tp, uj, and /2 are functions of r; as r — > co, we have 
to, p,, v — > and ip — > lnr. Note that this formulation is only 
sufficient for eccentricity resonances; if we were to consider 
inclination resonances, we would have to include the 0(z 2 ) 
terms in the metric. 

The contravariant components of this metric are: 



4>t 



-oje 

-2i> 

— 2fi 



~2 -2v 

■ uj e 



and 



(2) 



Equation |T} has a residual gauge degree of freedom 
in the sense that we may freely reparameterize r — > f(?). 
We fix this by requiring e' y+ ^ + ' 1 — r. This choice is easily 
verified to be valid for the Schwarzschild coordinate system 
in the case of a nonrotating black hole, and for the Boyer- 
Lindquist coordinate system in the case of a rotating black 
hole. 



2.2 Particle trajectories 

We utilize the Hamiltonian formulation of the equations of 
motion for a particle. As is well-known, the action for a par- 
ticle of mass p is S — — \i J dr, where r is the proper time 
along the particle trajectory. For our purposes, the fastest 
route to the torque formula is not to use the covariant rep- 
resentation of the action parameterized by the affine pa- 
rameter, but rather to explicitly parameterize the particle's 
trajectory using the coordinate time t, which is always possi- 
ble outside the outer horizon. This method, which explicitly 
keeps only the 3 physical degrees of freedom, is best suited 
to a perturbation analysis. The formulation of the problem 
is standard - the La grangian in coordinate tim e is the basis 
of the exposition bv llnfeld fc Pleba riski (1963), an d Hamil- 
tonianization of the coordinate time is a standard tec hnique 
in post-Newtonian calculations (e.g. lOhta et al.lll97&t ) - but 
will be explicitly given here since we will need to refer to it 
repeatedly throughout the calculation. 

Defining it M to be the 4-velocity, i.e. the forward- 
directed tangent vector to the particle's trajectory with 
u a u a = — 1, we see that dr = dt/u 1 and hence the La- 
grangian is 

L=^ = -± (3) 
dt ut w 

The degrees of freedom of the particle are its spatial coordi- 
nates x z (t); we note that u depends on the spatial coordi- 
nates x % and time coordinate t, and on the 3 spatial velocities 
x 1 — dx l /dt. The conjugate momenta are 7r; = dLjdx 1 . 

Noting that x' — u 1 /it* , we see that varying g a pu a vP — 
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— 1 at fixed x a gives gapu^du^ = 0, or u a 8u a — 0. There- 
fore, 



= utSu + Ui8(u x') = (itt + UiX % )5u + u mSx 1 

This implies that 

du l u f Ui u f Ui 

dx 1 «t + UiX % ut + UiU % I 'u l ' 



(4) 



(5) 



recalling that by normalization of the 4-velocity UiU 1 = 
— 1 — utu* shows this to be equal to (u*) 2 «i. Therefore the 
conjugate momentum associated with Eq. is 



dL fi t 2 



(6) 



where pi are the spatial components of the covariant physical 
4-momentum. From now on, we will simply write pi and drop 
the 7Ti notation. 

The Hamiltonian is given by 



H 



■ PiX 



u u 

u l u 1 



-pt, 



(7) 



where again we used the normalization of the 4-velocity, 
UiU 1 = — 1 — utu f . Thus the Hamiltonian for the particle's 
motion is simply the energy seen by an observer moving 
orthogonally to the hypersurface of constant t. From a dy- 
namical perspective, the Hamiltonian should be thought of 
as dependingon t, x z , and pv, the formula for pt is the mass- 
shell relatiorjj (derived from normalization of u): 



H(t,x\p t ) 



g tt g i 0p i p ] - y?gi 



(8) 



2.3 Nearly circular, equatorial orbits 

We now consider the nearly circular orbits in the background 
space-time. We restrict ourselves to equatorial orbits with 

z — p z = 0. 



2. 3. 1 Form of the Hamiltonian 

A circular orbit is a solution for which r = or (equivalently) 
p r — 0. We will be considering nearly circular orbits, i.e. we 
will expand the Hamiltonian to order (pr) 2 . From Eqs. J2]) 
and JSJ, we find that in general 



H = Cbp^+ eV^ 2 +e- 2 ' I (p r ) 2 +e- 2 *( W> ) 2 . 



(9) 



We now consider linear perturbations around a reference 
circular orbit. To do this, we first expand to second order in 

Pr- 

H{p^r,p r ) = g(p ,r,O)+ (P :> 

^/fi 2 +e- 2 ^{p 4> ) 2 2 

+0[( Pr )% (10) 



where 



H(j>^, r, 0) = Q Plll + eV^ + e-2*(p )2 



(11) 



For a given value of p$ one can find the minimum of 
H(p,f,,r,0) with respect to r,which (since dH/dr — 0) cor- 
responds to a circular orbit|j We can expand around any 
such minimum (with p$ = P$ and r = R) by writing 



Aptf, = p<p — P<t> and Ar = r — R. 



(12) 



The transformation from (r, (j),p r ,p$) to (Ar, 4>,p r , Ap^) is a 
simple translation and hence is canonical. Then H (p$, r,p r ) 
can then be Taylor-expanded around (P^,i?,0): 



where Cg ± 2 p 3 are the expansion coefficients. In order to 
study small perturbations of the orbits, we need to keep 
terms up to second order, i.e. /?i +/J2 +/?3 ^ 2, and we drop 
those whose coefficients vanish. This leaves us with 



H(Ap ll> , Ar,p r 



1 2 
Cooo + CiooAp^ + -C2oo(Ap^,) 

1 2 1 2 

+ -C 2oAr + -Co02(?5r) 



+CiioAp^AT\ 



(14) 



2.3.2 Relation of the coefficients to specific energy and 
angular momentum 

Some of the Taylor expansion coefficients in Eq. (| 14p have 
a straightforward interpretation, and all are calculable in 
terms of metric coefficients and the specific energy and an- 
gular momentum. We denote the specific energy (H/fi) and 
specific angular momentum (p<j,/fJ,) associated with a circu- 
lar orbit of radius r by £(r) and C(r). We may also define 
w to be the 4-velocity associated with the circular orbit. Its 
covariant components are wt = —£(R), w$ = C(R), and 
w r — w z = 0. Using the inverse metric, the contravariant 
components are 



tu' = e~ 2v {£ — ui£), 

= die" 2 " (£ — CjC) + e~ 2 ^£, 

and w r — w z = 0. 
By definition, 



Cooo = p£{R). 



(15) 
(16) 



(17) 



If we consider a sequence of circular orbits parameterized by 
r, we may take the total derivatives of the Hamiltonian with 
respect to r [i.e. derivatives in which p$ = /i£(r) varies as 
we take the derivative]: 



d 8H + ^8H 

dr or ar op<t> 

which since dH/dr = simplifies to 

£'{r) = C'{ r )^L 
0P4, 



Ci 



C'(R) 



n{R). 



(18) 



(19) 



(20) 



3 The requirement that the particle travel forward in time implies 
that we use the negative branch of the square root. 



4 We consider only the stable solutions; maxima of H or values of 
Ptf, for which there is no circular orbit solution are not of interest 
here. 
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We note that for a circular orbit, <f> — dH/dp^ = Cioo = 
Q(R), so Q,(R) can be interpreted as the angular frequency 
of the orbit as seen by a distant observer. For the circular 
orbit, we also see that trivially 



(21) 



Taking yet another total derivative of Eq. (|19p gives 



I d{r) f d 2 H j d[ M £(r)] d 2 H 
\ drdp^ dr <9(p^) 2 



(22) 



or at r — R: 

£"[R) = £"(R)!^± + c'(R) [Cno + fMC2 00 jC'(R)]. (23) 

Using the quotient rule, this can be expressed in terms of 
Q{R): 



Q'(R) = C 110 + nC2oo£'{R)- 



(24) 



and 



C020 = -fitf(R)£'(R) +iie- 2 "- 2i '{w t )- 3 C' 2 {R) 



(35) 



Explicit evaluation of these expressions is aided by a 



relation for w . Using the normalization g a pw a w 
and w — S7u>', we find 



P 



2v 2i> 



~,\2l-l/2 



(36) 



This completes the description of the C coefficients in 
terms of the commonly tabulated functions £(R), jC(R), and 
n(R). It is convenient also to define the epicyclic frequency 



k,(R) = V C020C0 



(37) 



it is easy to see that if Ap0 = 0, Eq. (|14[l guarantees that 
k(R) is the frequency of radial oscillations as measured by 
the coordinate time t. We also define the specific epicyclic 
impedance 



fiZ(R) 



C020 



Coo 



(38) 



Similarly, taking the total derivative of the relation 
dH /dr — gives 



8 2 H | d[/i£(r)] d 2 H _ Q 



(25) 



dr 2 dr drdp^ 

Evaluated at R, this simplifies to 

C 2o+fiJC'(R)C 110 =0. (26) 
A relation for C002 can be obtained from Eq. (|10[) : 



C002 



(27) 



The value of the square root is obtainable from Eq. 
giving C002 = e l/_2A1 /[e""(J/" - ojp^)], or 



C002 = 



£(R)-CjC(R) 



-2/3 



(28) 



Finally, we note that C200 can be obtained by directly 
taking the second partial derivative of Eq. noting that 

w, f, and do not depend on p^, we find 



9 2 # 



/iV" 2 * 



d{p^ 2 pr=0 
or at r — R: 



[ M 2 +e -2V( P(/> )2]3/2 



C200 = 



M -v- 2 * 



(29) 



(30) 



[i + e- 2 *r 2 (i?)] 3 / 2 ' 

Further simplifcation is possible if we apply Eq. (|11[> to a 
circular orbit, yielding 



£ = ii£ + e'Vl + e- 2,/ '£ 2 ; 

since £ — u)£ — e 2 "m', we conclude that 

-1 1 -2i/> ,-2 2v / t\2 

1 + e v £ — e (id ) . 

Substituting these results into Eq. (|30[) gives 

n —1 —2v—2tb/ t\— 3 

C200 = M e (^ ) 

Combining with Eqs. (|24[) and (|26|) gives 

dm = n'(fl) - e- 2 "- 2 *^*)- 3 ^^). 



(31) 
(32) 
(33) 
(34) 



3 RELATIVISTIC RESONANT TORQUE 
FORMULA: FORMAL SOLUTION 

3.1 Perturbation Hamiltonian 

Our next concern is the canonical treatment of a perturb- 
ing body. We separate the perturbation Hamiltonian into a 
perturbed and an unperturbed piece: 



H(t,x\pi) = H (t,x\pi) + Hi(t,x\pi 



(39) 



In Newtonian theory, the perturbing Hamiltonian Hi is sim- 
ply the gravitational potential of the perturbing body (plus 
an "indirect term" in formulations that do not use an iner- 
tial reference frame). In GR, there is a perturbation to the 
metric: 



(0) 1 u c>P 

9<*0 = 9 a p +h aB or g ' = g 
and Hi is then the variation of Eq. at fixed pi, 



,(0)a/3 



K 



>/3 



#1 = h° 



dg 



(40) 



(41) 



The latter can be obtained by varying the mass-shell rela- 
tion, g al3 p a p l 3 = -^ 2 : 

5g afi p a pp + 2g afi p a &pp 



■ 0. 



(42) 



Since Eq. ()41[) is defined at fixed Pi, the last term may be 
restricted to /3 = t, and: 

h" p p a pp 



6pt 



1 



og p a p/3 = 



(43) 



2p t " ^ p 2p t ' 

where we have used the rule that the variation of the con- 
travariant metric is —g aK g l3x 8g K \, i.e. the negative of h with 
its indices raised. Thus Eq. (141[l simplifies to 

h afi p a pp 



Ht = 



2p< 



(44) 



When doing perturbation theory, it is most convenient 
to do the explicit 3+1 expansion of the numerator and recall 
that 



t tt , ti 

P =9 Pt + 9 Pi 



tt TT I tl 

-g H + g pi 



(45) 
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H 1 = - 



h u H 2 + 2h u Hpi + h ij p iPj 
2(g u H ~ g^pi) ' 



(46) 



In first-order perturbation theory, it is permissible to replace 
H with Ho on the right-hand side of Eq. (I46p since the latter 
already has one explicit power of h. 

In the unperturbed case, the angular momentum p^, and 
the energy H are conserved. In the perturbed case, these 
variables change in accordance with 



P<t> = {P-P,H} P = {p$,Hi} P 



dHi 

' d<j> 



(47) 



and 



H = 



dH 



(48) 



_ dHi 

dt dt ' 

where {,}p represents a Poisson bracket. Since the pertur- 
bation arises from a secondary on a circular equatorial orbit, 
then the perturbation rotates at a pattern speed Q s given 
by the orbit of the secondary hole, i.e. Hi depends not on 
t and <j) individually but only on the combination (f> — Q, s t. 
This implies that the partial derivatives in Eqs. (|47|l and 
(|48[) differ by a factor of — f2 s , so 

(49) 



h = -a 



3.2 Effect of perturbation on the disk 

We consider an ensemble of particles initially in circular 
orbit at radius R with longitudes <f> equally distributed in 
4> £ [0, 2n). The perturbation Hamiltonian is assumed to 
turn on at time ti, and we wish to measure the torque 
on the disk of particles at some later time t?. The inter- 
val ti — t\ should be long compared with the orbital time 
f2 _1 , but short compared with the libration time so that 
first-order perturbation theory for the positions of the par- 
ticles is valid. It is apparent that the torque (T) averaged 
over the ensemble of particles must be second-order in h be- 
cause if Eq. (|47[) is averaged over the unperturbed particle 
trajectories, we find 



1 f 27r riH 



(50) 



In order to get a nonzero torque, we must compute the par- 
ticle positions to first order in perturbation theory, and then 
apply Eq. (|47[) . This will lead to a result that is second order 
in h. 

In what follows, we will construct the Green's function 
solution for the perturbations to the disk. In order to evalu- 
ate the late-time torque, we will decompose the perturbation 
Hamiltonian into Fourier modes in the longitude direction, 



Hi 



(m) 



m— — oo 



where each mode has the dependence 



(51) 



(52) 



on longitude and time (the latter is required since the per- 
turbation rotates with the orbit of the perturber). Since the 



Hamiltonian is real, H\ 



(m)* 



(-m) 



The torque transfer 



from different values of \m\ can be considered separately. 



This is because a first-order perturbation introduced by the 
m component will have an e"™* 1 '' 1 ' longitude dependence, 
and hence can only produce an azimuthally averaged torque 
when acted on by the m'th Fourier mode of the perturbation 
if m + m' = or m = —m'. 

Our final step in determining the first-order perturba- 
tion to the disk will be to integrate the Green's function over 
time, keeping only the resonant terms. 

3.2.1 Green's function solution for the perturbed disk 

We can compute the final position of a particle initially at 
4>i = <f>(ti) via a Green's function method. We consider first 
the effect of a 5- function perturbation at time t', i.e. we 
apply the perturbation 



W(t) = Hi(t)5(t-t') 



(53) 



Then to first order in perturbation theory, the perturbations 
to all variables can be written as an integral of the pertur- 
bation to that variable due to W over the range ti < t' < £2- 
Immediately prior to the application of W, the particle is at 
position 



<t>{t' - e) = <f>i + Sl(R) (t'-h). 



(54) 



Immediately after the application of W, any phase-space 
coordinate X undergoes a jump: 

OX = X(t' + e) - X(t' - e) = {X, Hi(t')} P . (55) 

One key difference between this and Newtonian perturbation 
theory is that since Hi depends on the momenta as well as 
the positions, the particle position can also undergo a jump. 
These jumps are: 

dHi(t') 



Or 

<>4> 

OPr 
Of0 



dp r 

dHijt') 
dp,^ 

dHijt') 

dr 
dHi(t') 



and 



(56) 



We then desire the final values of the positions and mo- 
menta. These can be freely propagated from t' + e using the 
unperturbed Hamiltonian, Eq. (|14J1 . The angular momen- 
tum is the easiest since it is conserved 
dHi(t') 



Apg>(t 2 ) = 0?V> = 



(57) 



The radial degrees of freedom are more subtle. They can be 
described by the Hamiltonian 

2 



H = iCo 20 (Ar+glH 
2 V 0020 



1 2 

+ -Cbo2(pr) + const, (58) 



which is identical to the Hamiltonian of a simple harmonic 
oscillator of effective spring constant C020, effective mass 
I/C002, and equilibrium position — (Ciio/Co2o)OfV- Under 
this Hamiltonian, the complex amplitude 



Ar 



Cno 
C020 



- , C002 

+ ll/7^ Pr 

O020 



(59) 



satisfies the equation of motion Z = — ik(R)Z and hence has 



-Ik(K) t 



dependence. Its initial value is 
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z(t' + e) - dHl ^ °} w dHl ^'- 



C002 dH x (t') 



C020 dr 



(60) 



dp r C020 d(p 

We may find Z(t 2 ) by multiplying by e~ iK(R) te-O. Taking 
the real and imaginary parts gives 

Cm, dH^t') 



Ar(t 2 ) = 



C020 d(p 

dHijt') Cho agi(t') 



+ 



\k(R) (h - 1')] 



dp r C020 0(f> 



and 

Ap r (t 2 



dH x (i?) 
dr 



MR) (fa - 1')] 



ggiftO _ Cno agiftp 

<9p r 



C) 



: /pH£ sinKi?) (fc -«')]. 

V (-'002 



(62) 



Finally, we may find the change in the longitude. The change 
in its rate of advance can be found by varying dHo / dp^, using 
Eq. (O, 



A0 = C2ooAp + CnoAr; 

and the change at time t 2 is found from 



(63) 



A<t>(t 2 ) = <><!> + A0(t 3 )dt 3 . (64) 

J t'+e 

Using Eqs. (|57[) and (|62|) . we may perform the integral to 
get 



A0(t 2 ) = 



C 2 oo(t2 — t') 7^7 — - 



, r [Cno dHig) , 



+ 



dHitf) C 110 dHi(t' 



8p r 



C020 d(p 



sin[«(J?) (t 2 - t')] 



k(R) 



C002 1 - cos[k(R) (t 2 - t')} 



C020 dr 



k(R) 



(65) 



3.2.2 Perturbed particle position for a particular Fourier 
mode of the perturbation 

At this point, we assume a particular Fourier mode m. 
Then the perturbation Hamiltonians have a dependence 
e im(^-Q a t). Since on the unperturbed trajectory, <j> advances 
at a rate Q(R), we may write 

dH^t 1 ) _ c im[si(R)-n B ](t'-t 2 ) dHi_(t 2 ) 
dr dr ' 

where the latter derivative is understood to be evaluated at 
the unperturbed longitude 

4 0) =& + n(.R)(ta-ti). (67) 

[Note that the actual longitude is 4>(t 2 ) — 2 ' + A(j)(t 2 ).] 
Inserting this dependence into Eqs. (I57p . (|61|) . (|62|) . and (|65p 
gives the following results for the J-function perturbation. 
For the angular momentum, 



(66) 



Ap<f,(t 2 ) = Op^ = — e 



dHi(t 2 ) 



(68) 



For the radial position, 

Ar(t 2 ) = ClW c im [ n ( fi )- n °Kt'-t2) ggiftg) 
C020 90 



+ 



affifcO Cno agiftg) 



9f>r Co20 9^ 

xe i m[ Q(fl)-n 3 ]( t '-t 2) cog[K(i?) (fe _ ^ 



C002 dHx(t 2 ) 



C020 dr 
xe i»[(i(jj)-n,](t'- {j ) s i n [«( J R) (t 2 - t')]. 

For the radial momentum, 



(69) 



A Pr (t 2 ) = _e-[«(«)-^]C*'-'=) C08[/s(iZ) (ta _ 



C020 
C002 



dHi(t') Cno SfTi(t a 



C020 



xe im[n(«)_n.](«'-i 9 ) sin[fi . (i?) (t2 _ 
For the longitude, 

A0(t 2 ) = 

+C110- 



(70) 



+C110 



Cup dHi(t 2 ) _ 
C020 90 



dp r Co20 

sin[ K (i?) (t 2 - *')] 



S (B) 



-Ci 



C002 9ifi(i 2 ) 1 - cos[ K (i?) (fa - f )] 



C020 <9r 
m[n(ij)-n B ](t'-t 2 ) 



(71) 



3. ,2. 3 Integration of resonant terms 

We now integrate Eqs. ()68H7ip over At' . There are many 
terms, however most are of short period. We therefore eval- 
uate only the Lindblad resonant terms, i.e. those that satisfy 
the condition 



m[Q,{R) - Q s ] « ±k(R). 



(72) 



For positive rn, the + sign is appropriate for interior res- 
onances and the — sign for exterior; for negative m this 
is reversed[f| It is convenient to write a resonant de-tuning 
function, 



D{R) = m[Q(R) - fi s ] =F K(R). 



(73) 



Within this resonance condition, we may replace the time 
integral involving cos[/t(ii) (t 2 — t')]: 



■ e i m [n(H)-n.](t'-t a ) cos[k{r) {t2 _ t i )]dt > 
2iD(R) ' 



(74) 



5 There are also corotation resonances where Q(R) f2 s , but we 
will not examine them here as the secondary hole actually orbits 
within the corotation resonance. 
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where At = t2—t\. A similar simplification occurs for the 
sine integral, 



,i m[ fi(fl)-n s ](t'-i 2) gin[fi . (jR) (<2 _ f ^ df > 



1 - e 



-iD(R)At 



2D(R) 



(75) 



Only these factors have resonant denominators. Integrating 
Eqs. (I68H71I) gives no change in the angular momentum, 

(76) 



A W (fa) = 0; 
for the radial displacement, 



Ar(t 2 ) = 



i-e- iD(H)A ' r .dH[ m) (t 2 ) 



+i 



2D(R) 

Cup dH[ m) (t 2 ) 
C020 dcj> 



C020 d; 
for the radial momentum, 

1 _ e -iO(fl)At 



(77) 



Ap r (t 2 ) = 



.affj m) (t 3 ) 



2D(i?) 



Co. 



dio 9#i (m) (t 2 



V C020C002 d(f> 
and for the longitude, 



(78) 



A0(t 2 ) = 



1 „-iD(K)At ^ 
1 — e Quo 

2D(i?) «(#) 

„ MT^ta) ^Ciio9^ m) (t 2 ) 
TO110 - ± 



C002 agf°(tg) 



Co 



dr 



dp r 



Co20 90 



(79) 



3.3 Net torque on the disk 



We are now ready to compute the torque exerted on the 
disk. Since the torque on the unperturbed disk vanishes, 
we may compute the (f>^ -averaged torque on the first-order 
perturbed disk. Recalling that only the — m component of 
the perturbation gives an angle-averaged torque on the m 
component of the perturbation, we find 



T 



(m) 



E 

» E 



d 2 H[- m \t 2 

dcj)dX 



■AX(t 2 



dX 



AX(t 2 



(80) 



X e{r,<£,Pr.P4,} 

Using Eqs. (|76H79[) . we may evaluate this as: 



T 



(m) 



l _ e - i D(R)At f d jj( 



2D(R) 



r (m)* r 



dr 



.dip dH { r\t 2 ) 

+1— hr- 1 — - ± 



. dH[ m \t 2 

8p r 

C^dH[ m \t 2 )' 



+ 



C020 90 V C020 dr 

0H[ m) * \.dH[ m) (t 2 ) [C^dH[ m) (t 2 ) 



8p r 



dr 



Co02 dp r 



C110 0H[ m) (t 2 



V C020C002 90 
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dH[ m) *(t 2 ) Cn 



+ 



90 k(R) 



C002 9#< m) (t 2 ) 
C020 9r 



dH^\t 2 ) 



C110 dH[ m \t 2 



C020 90 



(81) 



The quantity in braces {} looks complicated, but if we sub- 
it si 

(fa)l 



stitutc = VC020C002 (c.f. Eq. I3T[) it simplifies to 

1 _ e -iO(ii)At 



T (m) = ±im: 



2D(i?) V C002 
where the interaction amplitude is 



C020 i 5 ;m)/. 



(82) 



: (m) 



CllO rr (m) _ 



C020 



C002 
C020 9r 



(m) 



+ i 



95} 



(m) 



(m) 



1 a^i 



(m) 



/i£'(i?) ^(iJ) 9r 



+ i 



9i/J 



dp r 

(m) 



9Pr 



(83) 



We note that while 5 ( - m ' is formally evaluated at time t = t 2 , 
its time dependence is oc e _lmn =* anc [ hence its modulus 
jS (m) (t)| is constant. We can also relate 5 (m) to S { ~ m) : since 



H[ m ' = //} mJ * , and since the sign of the resonant term [c.f. 



r(m) 



Eq. J73}] changes when we switch from m to —m, 

S (-m) = _ <S (m)* > (g 4 ) 

We can then write the toiaZ torque arising from both the m 
and -m resonant terms as T = T (m) + T (_m) : 



T = ±i 



iD(fl)At 



2D(R) 



— (m o — m) 



C020 1 
C002 



,(m)|2 



(85) 



Since D(R) flips sign between the m and — m resonances, 



±i 



a -lD(il)At 



iD(R)At 



2D(R) 



+ ra- 



20 I <j(m) |2 



C002 



sin[£)(ie)At] /C020 



D(R) 



-2D(R) 



;(m)|2 



Coo; 



(86) 



Equation (|86[) has now separated into two pieces. There 
is an independent prefactor that contains the form of the 
resonance, and the factor 5^ m ' that encodes information on 
the normalization of the resonance and does not vary signif- 
icantly across its width. The first piece can be simplified by 
noting that it is dominated by regions with \D(R)\ < At~ . 
Its integral is 

sin[D(R) At] tt 

D(R) dR -W(Rj\> (87) 

so we approximate it near resonance as 

S m[D(R)At] tt 

~~ d(r) W(kJ\ S{R ~ Rl 

where R T is the radius of exact resonance. We thus find 
7rra „,^,„( m ) |2 



T 



\D'(R)\ 



/iZ(R)\S \ 5(R — R r 



(89) 



Often we want to know the torque density dT/dr. For 
a thin disk with proper surface density E, i.e. whose 3- 
dimensional proper density is po = £<5(z), the rest mass 
per unit radius is 



8 Hirata 



dr Ar J 



pow ■ nd V, 



(90) 



where d 3 V represents the volume of a spacelike 3-surface 
spanning the range from r to r + Ar and n is the unit 
forward-directed normal to this surface. Taking the surface 
to be at constant t, the normal is n a = (— e", 0, 0, 0) and the 
volume element is d 3 V = e* + ' i dr d<f>dz. This leads to the 
result 



^ = 2tiWE. 
dr 

It follows that 



dT 
dr 



(91) 



(92) 



4 RELATIVISTIC RESONANT TORQUE 
FORMULA: EVALUATION 

Having the formal solution for the torque (Eq. I89[l is only 
part of the problem; we also need the resonant amplitude 
g(m) Tjjjg section evaluates the amplitude and then shows 
that (within some restrictions) it is gauge-invariant. 



4.1 Evaluation of 5 (m) 

Here we require both the perturbation Hamiltonian, and its 
derivatives with respect to r and p r . These are all to be 
evaluated at the unperturbed circular orbit using Eq. (|46|) . 
For Hi itself, we see that since Ho ~ )J.£(R) and = 



fx 2 „ h u £ 2 - 2h t *£C + h^C 2 

Hl= 2° T^c • 



(93) 



For the partial derivatives with respect to the coordinates, 
we find that in general 



dHx 

dx k 



h u , k H 2 + 2h t \ k H oPl + /r 
2(g u H - g u pi) 

2h u Ho(dH /dx k ) + 2h u (dH /dx k ) Pi 

2(g*H-g* i p i ) 
„ g u (dH /dx k ) + g t \ k H - g l \ kPl 

-Hi 



(94) 



g u H - g u pi 

where the last term is associated with the derivative of the 
denominator in Eq. 1|46[) . Since dHo/dr — on a circular 
orbit, this implies 



dHi m 2v ti\ T £ 2 

dr ~ 2 6 



2fe^, r £:£ + fc^, r £ 2 
£ — CjC 



(95) 



For the partial derivatives with respect to the momenta, 
the general expression is 

dHi = (h u H + h u Pl )(dH /d Pk ) + h tk H + h ik Pi 
dp k g u H - g u pi 



Hi 



g u (dH /dp k ) 



(96) 



g tt: Ho — g u pi 

For the specific case of p r , we note that at the circular orbit 
dHo/dpr = and g tr = 0, so 



dHi 

d P r 



,{m) 



£ — uj£. 

We may now assemble the pieces to compute ,S^ m ' 

e 2u ( r m 1_ / UJ, r C 
2(£ -uC)\[ C T Z V V ' r + £ - CjC 



(97) 



41 



£ z - 2h 



{m)t(p 



£C + h 



1 



(m) " r £ 2 - 2h 



(m)tcj) 



, r £C + /i (m) ^, r £ 2 ] 



-ih 



(m)tr 



£ + ih 



(m)r<p 



(98) 



Note that this is independent of the particle mass \i and 
linear in the perturbation h al3 . 

It is possible to rewrite Eq. (|98[) in terms of the circular 
4- velocity w. Multiplying through by it)* = e~ 2 " '(£ — CjC) 
gives 



* c ( m ) 



w S 



m 1 , 

c : l [ 

" 2Z 



£ — CjC 



+ ih 



(m)ro 



h^ m)af3 w a w 



(99) 



This form will be most useful in proving gauge invariance 
and in practical applications. 

4.2 Gauge invariance 

In general the perturbation h a/3 could be expressed in many 
choices of gauge. These differ by the relation 

ha/3 — > h a p — ^a-p — £j3;ae (100) 

Since Eq. (|99[> is linear in h a @, the contributions to 5' m ' 



from the pre-existing and gauge perturbations simply add, 
so to show invariance of the torque it is sufficient to prove 
that a pure gauge perturbation 

h a fi = — Ca;/3 ~ £/3;a (101) 

leads to zero resonant amplitude S^K We restrict our at- 
tention to gauges that preserve the fundamental symmetries 
of the problem, i.e. that have reflection across the equatorial 
plane and have helical symmetry, where the m Fourier com- 
ponent has an oscillatory time dependence oc e _in = mt . With- 
out loss of generality, we may consider the Fourier modes one 
at a time, so we will consider the order m Fourier mode be- 
low and avoid writing the superscript ■ rn ' explicitly. Further- 
more, it is easily seen that the z coordinate is superfluous 
in computing Eq. (|99p in the equatorial plane, so we may 
restrict ourselves to the 2+1 dimensional equatorial slice of 
the spacetime. 

While one could solve for for a pure gauge mode by 
explicit evaluation of Eq. (|1U1[I followed by substitution into 
Eq. (|99[) . it is far easier to solve the problem by defining the 
combinations of metric perturbations and 4-velocities that 
appear in Eq. (l9"9l : 



h 
Is 



i a B 
-n a f3W MT , 



— h a ^ >r w a Wf3, and 



h ra w a , (102) 

and evaluating these in terms of £ with the help of Lie deriva- 
tives. We may then substitute into 
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w S — 



m 1 / W, r £ 

=F — l 2v r H ! 



hT^+ih, (103) 



and then check whether the terms add to zero. 



4-2.1 Evaluation of h 

We begin by writing the equation for h a p (Eq. I101[) in the 
al ternative fo rm using the Lie derivative (e.g. Appendix C 
of lWaldll983, 



h a /3 = — £^g a a, 

h af) = £^ = ?g^ n - g^e,y ~ S°V,7- 



(104) 



(105) 



We are now in a position to compute the required term 
hapw^vjP . Noting that g a pw a w 13 = —1, we take the Lie 
derivative £$, 

— h ct pw a w 13 + g a f){££w a )w f> + g a 0W a (£^w 13 ) = 0; (106) 
using the symmetry in a and j3 then gives 

^h al 3W a vj 13 = g af jW°'(£ 6 vj 13 ) = w{£^w) =w-[£, w], (107) 

where [,] denotes a vector commutator. We explicitly eval- 
uate the ' and * components of the commutator; recalling 
that w r = and iu* and depend only on r, we find 



rj. it s-r t t J-t <b t~t 



(108) 



using SI = (j> = wr/tu* and the angular and time depen- 
dences of £, we conclude that 

[$, wf = f«/, r + imio l (fl s - fi)£*. (109) 
Similarly, 

[$, wf = e r i»*r + imw*(O s - fi)^. (110) 
Taking the dot product with w gives 

W-[£,w] = ^ r (w t W t ^r + W^W.r) 

+im(O s - Q)™*^' + tw^*). (HI) 
The first term evaluates to zero: 



«;*(£' - OC') = 0. 



(112) 



(The first equality can be shown by differentiating the rela- 
tion w a w a = —1 with respect to r.) The second is simplified 
using wt = — £ and Wcf, = jC; thus we find that in general, 

h = i/i a/3 w Q u>' 3 = im(a - + ££ ). (113) 



6 The Lie derivative does not generally allow raising or lower- 
ing indices. The raised-index relation arises by considering the 
inverse-metric formula g a ^gp 1 = 5%. Applying the product rule 
gives (£^g aP )gp 1 + g a P£ i gp 1 = 0. Substituting Eq. iTTOi} and 
raising indices then gives £^g a ^ = h a P . 



4-2.2 Evaluation of I2 

We now turn our attention to I2, which appears in the sec- 
ond term in Eq. (|f U3p . A reorganization gives 



h = ( ^h al3 w a wp ) - h afi w a Wf). 



(114) 



Use of Eq. (fTOl) gives 

^—h a ^'w a u!0 S j = im[— Sl'w* + (Sl s — Sl)w t t r] 

x(-£? + £(+) 
+im{Sl B — Sl)w 

x (-£'? + C'^ - £Z\r + C^, r ).(115) 

To complete the evaluation of I2, we introduce the 1- 
form field 



S/3 = £g/g r W/3, 



(116) 



whose components are sp = lup >T , or explicitly s t = —£', 
= £' , and s r = 0. Then the last term in Eq. (|f f4p is 
—h al3 WaSp. We can see that 

g al3 w a Sf3 = w 13 ' sp — —w t £' + C! = (H7) 

since SI = £'/£' = w' l '/w t . Taking the Lie derivative £ s 
gives 

= h al3 w a sp + s a £^w a + w a £^s a . (118) 
Rearranging and expanding the Lie derivatives gives 

-h ^WaSp = S ^W ai p+S Wp^, a +W f Sc,/3 



1 a j-P 



(119) 



Recalling that the terms containing w a: /3 and s a ,p are only 
nonzero for f3 = r, that = fill?*, and the oc e im(0 - nst) 
dependence of the components of £, we reduce this to 



h w a sp 



+im(s' 1 ' - n s s t )(-£Z t + 



+im(n - fi s )to*(-£'£' + C'^). (120) 
Combining this with Eq. (|115|l gives 
h = im[-Q.'w t + {Sis - Q)w\ r + - fi s s'] 

x(-£:c t + rr) 

+im(fi s - fi)w*(-^',r + ££*,r) 

+ - w'f" + ^ £")C r - (121) 

Further simplification of this equation is possible using 
the contravariant components of s: raising indices gives 

s = e {t — loL ) and 

= Cue- 2 " {£' - CjC') + e - 2V, £'. (122) 

From this we obtain 

-s'E' + s^C' = -e- 2u {£'-cj£) 2 +e- 2 *C' 2 
= £' 2 [-e- 2 "(n-w) 2 +e- 2 *] 
= e- 2y -^(™')- 2 £' 2 , (123) 

where the second line used SI — £'/C' and the third line 
used Eq. (|36[) . However, we also see that: 

t- W*Co20 

ui Zk = 
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-w il L + e (id ) L 



-w'Q'C' - s£' + s^C' 
t ££"-£'£' 



s£' + s^C 



<t> r' 



-w t + w L — sc + s vL 



(124) 



[Here the second line used Eq. (|35|l ; the third line used 
Eq. (I123|l ; the fourth line used Q = £'/£.' and the quotient 
rule; and the fifth line used that u>^/«/ = fl — £'/£'.] Equa- 
tion (I124|l leads to two major simplifications in Eq. (|12ip . 
The term involving £ r simplifies dramatically. Also, using 
the first and third lines of Eq. (I124|l and £' = Q.C' , we find 
that 

- nV + - ^s* = (n - n s )s' + w/;'" 1 ^^ (125) 

Therefore, Eq. (|121|l simplifies to 

h = im[(n s - Q.)(w\ r ~ s t )+w t £'- 1 ZK](-££ i + C^) 

+im(n s - n)w\-£^, r + c^, r ) 

WZkC- (126) 

A final level of simplification involves w*, r — s*. Using 
the explicit expressions, Eq. Q15p for w and Eq. (|122p for 
s', we see that 



s = e 



[— 2v >r (£ — CjC) + £ — Cj^C — U>C\ 

"{£' -QC 1 ). 

Cj r C 



-\2Vr- 



£ — CjC 

This allows us to eliminate s from our expression for I2 
I2 = im 



(127) 



-(n s - n) 2^ 



+w t C' Zk 



£ — CjC 



+im(a - ft)w*(-££* ir + £^ 



(128) 



^.l?.<? Evaluation of I3 

Finally, we consider ^3 = h ra w a . This is most easily com- 
puted by explicit evaluation of the contravariant components 
using Eq. (|1U5I) : 

h rt = -e- 2A £*,r + im(Cj - O s )e _2 T (129) 
and 

h T * = -e" 2/i ^,r + im(w - a)e- 2l/ cDC r - ime" 2 ^ r . (130) 
This implies 

h ra w a = e~ 2ii (£ £* >r — £^ ,r) 

+im(u — fi s )e~ 2l/ £ r (— £ + CjC) 

-ime" 2 *^''. (131) 



The terms involving £ r can be simplified using Eqs. ()15[) 
and (|16l) . which simplifies them to im^sw' — w^)£ r . Further 
using w 1 ^ = fitu* gives 

h = h ra w Q = e~ 2p -(£t,r - C^, r ) + im(ft s - ■ (132) 

The other contributions to S do not explicitly contain fi, 
so in order to prove gauge invariance we will need to elimi- 
nate fi in favour of other variables. Equation ((28)) provides 



a convenient choice: it and the definitions of k and Z tell us 
that 



— 2fi t r , t K 

e ^ = fiw G002 = w — . 

We thus arrive at our final expression for 73: 



(133) 



(134) 



^.^.^ Putting it all together 

We now substitute Ji, I2, and J3 into Eq. (|103[) . giving 
1 



w S 



m 1 / u> r C 

■a^z\ 2v - + £±QC 

xw\-£^ +C^) 



im(fi s — Q) 



. m 

T1 z 



-(n s - ft) 2//, r - 



Cj. r C 

£ — CjC 



+ - 



w 1 Zk 



c 



{-£C + cc 



T i-(a - n)to*(-£^, P + ce, r ) t w*k? 



(135) 



We may divide through by w f on both sides, and cancel 
the terms involving 2v, r — Cu, r C/{£ — CjC). Collecting the 
remaining terms gives 



5 = i£{-£e+£?)[m(n- 



+[m(n ~ n s ) =f • 



(136) 



In general, this is nonzero. However, there is one piece of 
information we have not used: that the resonant amplitude 
is to be evaluated at the resonance location D(R) = 0, i.e. 



m(fi — fi a ) = ±k. 



(137) 



When - and only when - we use this fact, we see that 
Eq. (|136[) vanishes. That is, the resonant amplitude 5' m ' 
is only gauge-invariant when evaluated at the resonant posi- 
tion^. This is not a problem since the torque formula contains 
a 5-function at the resonance. 

Thus we see that a pure gauge perturbation leads to 
zero contribution to S' m ' at resonance, and the resonant 
torque is gauge-invariant. 



4.3 Epicyclic geodesic formulation 

There is an alternative way of writing the resonant ampli- 
tude that will be better suited to computation in the 
Schwarzschild and Kerr spacetimes. We will argue in this 
section that is related to a particular integral of the 
metric perturbation along the world line of a test particle 
on an orbit with very small eccentricity. This formulation 
has some utility in the Newtonian case, but it will be shown 
to be very powerful in Paper II, where we will relate it to the 
gravitational waveform emitted by a test particle on such an 
orbit. It will thus allow computation of 5' m ' using standard 
methods for computing waveforms, without explicit evalua- 
tion of the metric perturbations. 
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Our starting point is to consider a particle on an un- 
perturbed orbit (i.e. traveling according to Ho) oscillating 
between r — R—e and R+e. To first order in e, its trajectory 
is given bj0 

r = R + e cos(nt), 
p r = /j,eZsm(Kt), 



— Qt + g ^ 110 sin(fi;t), and 



(138) 



Now we consider the integral of the metric perturbation 
over the test particle world line 



/ 

J V 



(139) 



where the range of integration is over any epicyclic period, 
i.e. from ti < t < ti + 2-k/k, for any t\\ and in the second 
integral, T a/3 is the stress-energy tensor associated with the 
test particle and V is the region of 4-volume in this range of 
coordinate time. By construction, It is linear in the metric 
perturbation. It is also invariant under gauge transforma- 
tions respecting the helical symmetry ' oc e im ^ _n, *' ) 
since under a gauge transformation, Eq. ()139p changes by 



eAI T = 



M Jv 



^T^^del^d 4 *. 



(140) 



We may integrate by parts to move the l( g derivative onto 
r a/3 ; but T ,a/3 ;/ 3 = for a test particle traveling along a 
geodesic. The boundary terms at t — ti and t\ + 2n/n also 
cancel each other since both £i m) and T Q/3 are invariant 
under translation in time and longitude by t — > t + 2tv/k 
and 

2tt 2tt 2tt 

4> ^ 4> + Q — = 4> + Q B — =F — > (141) 
k k m 

respectively. Thus A/t = 0, and 7t is gauge-invariant. 

We may use the gauge invariance of Jt and <S' m ': if a 
relation between them can be demonstrated in one gauge, 
then it must be valid in any gauge. We choose the gauge with 
h u = h t4, = h?* = 0. This gauge exists for the generic case 
where m^O, since one may write the gauge transformation 
relations 



, (142) 



Ah" \ 




( g u ,r "2c 




( e 




H 




^ 




Ah** j 




V 9**,r 


2c t j 





where we define the vector c by c t = Qsg 1 * 
g t * — £l s g u , and c r = c z = 0. If this 3x3 matrix A is 
nonsingular (which may be easily verified for some cases 
such as Schwarzschild), then the gauge h u = h l * = h** = 
exists everywhere. (We will remove the condition on A later.) 
In this gauge, we find from Eq. Q99p 



7 These equations may be obtained from the Green's function 
relations in Section 13.2.11 by taking a particle in a circular orbit 
at radius R that passes longitude rfr = at t = 0, applying a 
perturbation at time t = that increments r by = e, and 
considering the solution at t > 0. 



5 (m) = i 



(143) 



We also find that in computing It only h rt and h r <j> con- 
tribute, and since u r is already 0(e) we find u r 
and 



-2£, 



It = 2c 



I smUt) — At. 



(144) 



We may raise r in the perturbation using the factor of e 2 ^ , 
and simplify this to 



-2\ZS 



Cm) | 



sin(^)e imW - nst) dt, 



(145) 



where the |o reminds us to evaluate S^ m ' at (f> — t — 0. On 
resonance the complex exponential decomposition of the sine 
allows us to evaluate the integral to —m/K, so we conclude 
that 

Z 



It 



-27r-S (m) . 



(146) 



So long as det this relation must be valid in all gauges 

since both sides are gauge-invariant. 

In some spacetimes there are radii where det A = 0; 
however we may show Eq. ()146[) to be valid there as well. 
We may consider a family of spacetimes M(P) whose met- 
ric tensor components are analytic in the parameter P, the 
desired spacetime is M(0), and detA oc P n (n — 1, 2, or 
3) for small P. Then we may carry through the argument 
for slightly different values of the parameters controlling the 
spacetime and prove Eq. (|146[1 ; then since both and It 
are analytic and equal in a neighborhood of P — they must 
be equal at P — 0. Thus Eq. (|146[) remains valid regardless 
of whether detA = or not. 

We may then express S' m ' in terms of the integral of 
the metric perturbation against the stress-energy tensor of 
a test particle on a slightly eccentric orbit, 



2^Zj/$ faB ^^ d4x - 



(147) 



A further simplification occurs if we extract the Fourier 
mode of frequency — mf2 s from the stress-energy tensor, 
which is the only one that can lead to a nonzero integral 
against h^ m ' . The t- integral is then trivial and we find 



: (m) 



efiZ 



% I h ( ^>f [ - mn ^ ^f^A+t~gA i x. (148) 
2 J 



A second version of the epicyclic formulation is as fol- 
lows. We note that the average amount of power V trans- 
ferred to the test particle by the perturbation in one radial 
cycle is 

r tt+2n/ K 

V= 2^J t Hl dt ' (149) 

the integrand can be evaluated along the unperturbed tra- 
jectory since Hi is already first-order, and we note that the 
dot pulls down a factor of — imQ s : 

*1+27t/k 

(150) 



V 



imfi s /iiK 



2tt 



2u l 



the integral simplifies to helT, and so we conclude that 



r = ; elT- 



(151) 
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From this, we extract a relation between and the power 
provided to a particle on a slightly eccentric orbit: 

2i 



mfL,/i.Ze 



V. 



(152) 



Note that V pertains to the particlar m-mode and hence 
may be complex if the peak power occurs at a resonant phase 
other than or 7T. 



5 NEWTONIAN KEPLERIAN LIMIT 

We now consider the limit of Eq. ()89p for non-relativistic 
Newtonian-Keplerian disks, and show that it reduces to the 
familiar result. 

In the limit of M/r <C 1 where we expect to recover the 
Newtonian result, the metric for a central object of mass 
M has v — -Mr" 1 < 1, ip = lnr, and u) = p, = 0. The 
Hamiltonian evaluated at zero radial momentum, Eq. (|lip . 
is 



- _ M + (p ) 2 - 



minimizing over r gives r 



2p?r 2 

(p^) 2 /(M M 2 ) = C 2 /M, so 



(153) 



C = M 



1/2 1/2 



(154) 



The energy is obtained by substituting back into Eq. (|153 
2r 

Using the results from Section [2.3.21 we find 



(155) 



n = M 1/2 r /2 . (156) 

Then w* = e~ 2 "{£-CbC) = l + §Mr _1 , so we find to leading 
order 

C002 = — and C020 = (157) 

This implies the epicyclic frequency and specific impedance 

«: = M 1/2 t-- 3/2 and Z = M 1/2 r~ 3/2 . (158) 

Since k = Q. in the Newtonian Keplerian case, the Lind- 
blad resonance condition m(Q. — Q B ) — ±/c is satisfied for 



m =p 1 



-Qs or r 



m =F 1 



2/3 



(159) 



We label the resonances with positive m so that the lower 
sign corresponds to the OLRs and the upper sign to the 
ILRs. In the Newtonian Keplerian case, there exist OLRs 
for each positive integer m, while the ILRs exist only for 
m ^ 2. 

The resonant torque further involves the detuning func- 
tion dD/dr: 

3 



D' = mO' =f «' = (m T 1)^' = 



(BIT l)M 1/2 r" 5/2 . (160) 



We now consider the resonant amplitude S' m '. In the 
nonrelativistic limit, the time-time metric coefficient is 



2$, 



(161) 



where $ is the Newtonian gravitational potential associated 
with the perturbation, and the other components are small0 
Keeping the leading-order (r 1//2 /i") terms in Eq. (|98p gives 



: (m) 



M- 1/2 [~2mr 1/2 ¥ m) T r 3/2 ^ m \ 



(162) 

The Keplerian analogue of the binary black hole case is 
for the perturber to be a point particle with mass qM, where 
the mass ratio q <C 1. Without loss of generality we may 
place the perturber at longitude A = 0; any other choice of 
longitude would result in $ ( - m ' and hence S^™ 1 ' being mul- 
tiplied by a factor of e~ lmA , which will have no effect on 
liS^" 1 - 1 ! 2 or on the torque formula. 

If we place this particle at radius r s , its perturbing po- 
tential is 



$(r, 



qM 



2r s r cos < 



H — j r cos < 



(163) 



where the second term is the "indirect" term resulting from 
the acceleration of the primary (i.e. it is necessary to keep 
the primary at the centre of the coordinate system, as is 
standard practice in celestial mechanics). We may project 
out the order m Fourier coefficient, 



$ (m) (r,0) 





qM 
2^7 



^e- im ^<J>(r,0) 
[-&$(?) + C <W], 



(164) 



where c = r/r s — [(m=p l)/ m] 2 ' 3 and b represents a L aplace 
coefficient (e.g. Eq. 6.67 of lMurrav fc DermotdkoOGv ). Then 
S' m ' (evaluated at zero longitude) is 



; (m) 



qM 



1/2 r 



2r ; 



1/2 



2m^ 2 b^l{,)~2e' 2 8 ml 



(165) 



The 5 m i term exists only for the m — 1 OLR (lower sign), 
so we may simplify this to 



_ {m) _ gM^y/ 2 

5 -~2l^~ 



S5 m l ± d>i?2 (?) 



(166) 



T 



Substitution into Eq. (|89|) then gives 

27T m q 2 Ma _ , 

T— — rM r— ; °\ r — -Rr) 



3 m =F 1 4r B 
x \2rnb[^^)-,5 ml ±^'^)\ 2 . 

Using m/(m =p 1) = ?~ 3,/2 , we reduce this to 
T = =F |w 2 M? 1/2 5(r - A) 



(167) 



(168) 



In terms of the disk surface density, E = fiS(r — R T ) /(2nr), 
this becomes 

7rVa 1/2 
T = =F^-F rME 



(169) 



8 For r/r s of order unity, £/£ ~ (rs/M) 1 ' 2 , but the components 
such as for nonrelativistic perturbers are suppressed by higher 
powers of M/r B . 
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At m ^> 1 or |<; - — 1| <C 1 we may meaningfully consider 
the smoothed torque density over many resonances(f] Noting 
that m(l — ? ) = ±|, w e find using the large m expa nsion of 
the Laplace coefficient (|Goldreich fc Tremaine|[l980r i: 

6^( f )«^o(|) and 6^'(0«±^i(|). (170) 
Then for large m, Eq. (|169[) becomes 

T = T^V^MS [2*0(1) + *i(|)] 2 • (171) 

This is for a single resonance. For a continuum of resonances, 
we need to substitute the resonance order m = 2r s / (3|r a — r\) 
and multiply by the density of resoances |dm/dr| to get the 
torque density 



dT 32 g 2 r E 4 MS 



2M 2 



[2*o(f) + *i(|)] 



(172) 



dr 81 (r s - r) 4 
which agrees with Eq. (18) of lGoldreich fc Tremainel (|l98(t ). 



6 RELATIVISTIC DISC HEATING AND 
SURFACE BRIGHTNESS 

Thus far, we have considered the angular momentum and 
energy transfer to the disk at the Lindblad resonances. In 
Newtonian thin-disk problems, it is often the case that the 
disc can radiate energy but not angular momentum. Since 
an orbit of fixed angular momentum has a minimum possi- 
ble energy, one can then compute the rate of energy input 
that does not go into orbital energy; this amount of energy 
goes into epicyclic motions, which are eventually converted 
to heat and ultimately radiated. The relativistic case is far 
more complicated because radiation carries away both en- 
ergy and angular momentum. We shall consider the problem 
here under the following two simplifying assumptions: 

• The dissipative process is localized, i.e. the energy of 
epicyclic motions is dissipated near the resonant radius 
rather than being transmitted to a distant part of the disc 
(e.g. via density waves). 

• The energy is radiated away locally, i.e. we assume a 
thin disk rather than an ADAF or other radiatively ineffi- 
cient solution. 

The second assumption is necessary in order to main- 
tain a thin disc, i.e. for the consistency of this paper. In some 
ca ses, it may w e ll bre ak down. For example, in the problem 
of IChang et ai] (|2010l ). in which the secondary "shepherds" 
the inner disc to smaller radii, it is conceivable that heating 
from resonant torques could destroy the thin disc solution. 
Even in this case, however, we would like to know the res- 
onant heating formula for a thin disc: inability to produce 
the required flux F for any disc temperature would be a 
sufficient condition for the destruction of the thin disc. 



6.1 Definitions and mathematical relations 



We use the formalism of I Page fc Thornel (|1974| ) to investi- 
gate the flux emerging from the disk, although we do not 

9 Whether dT/dr is really smooth depends on the nature of the 
dissipation mechanism, which we do not consider here. 



make the assumption that the disc is time-steady. We do 
assume that the disc is thin and that the internal energy 
is negligible compared to the orbital energy, i.e. if the bulk 
4- velocity of the baryonic material is u, 



rpfl 



pou u + v + u q 



q u , 



(173) 



where po is the rest mass density (i.e. the mass of a baryon 
times the number density), q is the heat flux, and t^ v is the 
stress tensor in the baryon rest frame (by definition /m u = 
and t M "it M = 0). In accordance with lPage fc Thornel l| 19741) . 
we assume that q lies in the z-direction (i.e. q z is the only 
nonzero component). The disc is assumed to be contained 
within a vertical thickness of \z\ < H; the stress tensor at 
z = ±H is assumed to satisfy 



t r — tt 



0. 



(174) 



iPage fc Thornel (| 19741 ) explicitly write time and longitude 
averages of these quantities, with the idea being to treat e.g. 
turbulent stresses as part of t ll v rather than as small-scale 
structure in u. We will not write these averages explicitly, 
but note that (i) they are implied, and (ii) in our case, the 
time averaging is assumed to be over a duration long com- 
pared with the turnover time of turbulent eddies but short 
compared to the evolution timescales of the system (e.g. the 
merger timescale). It is assumed that the disc material is 
on nearly circular orbits, but possibly with a small radial 
velocity, i.e. u l = w l , u = w , u z = 0, and \u r \ < y/M/r. 

We define the integrated quantities through the disc: 
the surface density, 



£(r,t) = / po(r,t,z)dz, 



and the integrated shear stress, 



W/(r,t) = / t/(r,t,z)dz. 



We also define the one-sided emergent flux 



F(r, t) = q z {r,t,z = H) = -q z {r,t,z 



-H), 



(175) 



(176) 



(177) 



which is the flux that would be seen by an observer sit- 
ting at the disc photo sphere and corotating with the disc 
l|Page fc Thornelll974l ). 

We further neglect stresses in the tangential direction, 
i.e. we set t^ = 0. Orthogonality with u then implies t^ — 
0. We note that the requirement that t^uF" — 0, combined 
with the approximation that u w r, gives us the integral 



t t r {r,t,z)dz = -QW^ r (r,t). (178) 



6.2 Conservation laws 

As is the case with the time-steady thin accretion disk, it 
is convenient to use the conservation of baryonic rest mass, 
angular momentum, and energy to solve for the state of the 
system. In our case, the equations will be time-dependent 
but their derivation is similar. For any current j satisfying 



J ;<x = r, 



(179) 
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where the source term F is the amount of charge added per 
unit proper 4-volume, we 

0=^(rj Q )+rr; (180) 
integrating this equation over z from —H to H, we get: 

= S- ( r [ f da j + ( r / f dz 



dr 



Of 



+r\j z (z = H)-f(z = -H)]+r Tdz. (181) 

J-H 

This implies, for the rest-mass current ' m 'j M = poii M , which 
has no source, 

= -^-(r'Eu) + rw t t. (182) 
Or 

For the angular momentum current < - L 'j M = , there is a 
source, namely there is an angular momentum dT/dr added 
per unit radial coordinate per unit coordinate time. This is 
related to the source via 

H 

= 2~r I (L) rdz. (183) 

H 



dT 



dr 
Thus 

= ^-(rWV + r£E H r ) + r£w t t + 2r£F + (184) 
Or 2k dr 

For the energy current ( - E - ) j M = — Tt M , the source differs from 
the angular momentum source in that the energy added is 
equal to Q a times the angular momentum added. This is a 
direct consequence of the fact that the time dependence of 
the metric perturbation consists solely of a pattern speed 
fi s . Then (E) T = n s (L) r, so 

= JLfrSlWf + r£T,u r ) + rEvft + 2r£F + — (185) 
Or 2-k dr 

Equations (|182l) . (|184l) . and (|182p provide 3 constraints 
for 4 unknowns (£, u T , F, and W^ r ). They can be solved if 
a prescription i s available for the shear s tress W ( p r , e.g. an 
a-prescription (|Shakura k, Sunvaevlll97 3')P1 



6.3 No-viscosity solution 

A special case of interest to us is the case where the viscosity 
of the disc is negligible (W^, r = 0). This limit is appropri- 
ate in the final stages of a binary black hole inspiral where 
the viscous timescale beco mes short compared t o the merger 
timescale, as occurs in the IChang et ail l^Old ) calculation. 
Then the disc evolution is dominated by angular momen- 
tum transport via the resonances and by the inspiral of the 
secondary black hole (itself driven by radiation reaction). 

Writing Eq. (| 184p without the W,f, r term, and using 
Eq. (|182[) to eliminate S, we find 



We have used the general expression for the divcrg 
( — 1//2 9a[(— Iffl) 1 / 2 ,?'™], and recalled that for our choice of co- 
or dinates — | <7 1 1 / 2 = r. 

11 I Page fc Thornd {l973) were able to solve this system in the 
time-steady case without assuming any prescription for angular 
momentum transport by setting £ — > and using the 3 equa- 
tions to solve for the remaining unknowns u r , F, and W^ r . This 
method is clearly not applicable to a transient event such as an 
inspiral. 



1 HT 

rC'T.u r + 2rCF = 

2-k dr 

Similarly, using Eq. (|185[) gives 

fi s dT 
2-7T dr 

This gives us a linear system for u T and F, with solution 
1 dT £ - ft s C 



r£'T,u + 2r£F = 



and 



F = 



27rrE dr ££' - ££' 



1 dT n s c - £' 



(186) 

(187) 
ion: 
(188) 



47rr dr ££' - ££' ' 



(189) 



To proceed further, we use Eq. (|92[) in the flux equation. 
We then reduce this using the relations 



££'-££' =£'{£-Q,£) 
and 

n s £' -£' = £'(Qs -Q) = T — , 

m 

yielding finally 

KK W l 



F = 



2r\D'\ £-fl£ 



Z-£\S {m) \ 2 5(r-Rr 



(190) 



(191) 



(192) 



Thus the emerging flux is, as expected, proportional to the 
surface density of material at resonance and localized at the 
resonance. In reality, the 5-function would be smeared out in 
a way that depends on the dissipation mechanism. We note 
further that r is not a proper radial coordinate: an observer 
sitting on the disc would measure a proper radial distance 
element e M dr instead of dr. That is, the emitted flux per 
unit length (units: ergs -1 cm" 1 ) along the circumference as 
measured locally by an observer on the disc would be 



F dr 



-KKe 11 2v „( m ) j2 



2r\D' 



-ZEIS 1 * 



(193) 



Equation (|193|l gives the emitted flux required for the 
disc to remain thin. It is of course emitted over some finite 
range of radii: there is a finite damping region for the density 
waves excited at each Lindblad resonance, and turbulent 
diffusion may transfer some heat to neighboring parts of the 
disc. If this amount of flux cannot be radiated by any viable 
disc model regardless of the temperature then the thin disc 
solution must fail. 



7 SUMMARY 

This paper has worked out the general formula for the torque 
on an equatorial disc in a stationary, axisymmetric space- 
time with an equatorial plane of symmetry due to Lindblad 
resonances associated with a perturbation. We have shown 
that the torque formula is gauge- invariant, and that the fa- 
miliar formula is recovered for the problem of a Newtonian 
Keplerian disc with a perturber on a circular orbit. We have 
also obtained the expression for the radiated flux required 
to maintain a thin disc solution. 

The most important astrophysical application of the rel- 
ativistic torque formula is to the Schwarzschild and Kerr 
spacetimes. The computation of the resonance locations and 
amplitudes for these cases is presented in the companion pa- 
per, Paper II. 
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